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We discuss the phase coherence properties of a mesoscopic normal ring 
coupled to an electric environment via Coulomb interactions. This system can 
be mapped onto the Caldeira-Leggett model with a flux dependent tunneling 
amplitude. We show that depending on the strength of the coupling between 
the ring and the environment the free energy can be obtained either by a Bethe 
ansatz approach (for weak coupling) or by using a perturbative expression 
which we derive here (in the case of strong coupling). We demonstrate that the 
zero-point fluctuations of the environment can strongly suppress the persistent 
current of the ring below its value in the absence of the environment. This 
is an indication that the equilibrium fluctuations in the environment disturb 
the coherence of the wave functions in the ring. Moreover the influence of 
quantum fluctuations can induce symmetry breaking seen in the polarization 
of the ring and in the flux induced capacitance. 

Key Words: persistent currents; Coulomb blockade; single-electron tunneling; 
quantum coherence 



1. INTRODUCTION 

Mesoscopic physics is concerned with systems that are so small that the 
wave nature of the electrons becomes manifest. The quantum mechanical 
coherence properties of small electrical systems are thus a central topic of 
this field. In this work, we are interested in the coherence properties of the 
ground state of a mesoscopic system, namely a ring subject to an Aharonov- 
Bohm flux, which is capacitively coupled to an external environment which 
exhibits zero-point polarization fluctuations [1] . In a mesoscopic ring, quan- 
tum coherent electron motion leads, in the presence of an Aharonov-Bohm 
flux, to a ground state that supports a persistent current [2-6] . The persis- 
tent current is thus a signature of the coherence properties of the ground 
state of a mesoscopic ring. Only electrons whose wave function reaches 
around the ring contribute to the persistent current. Thermal fluctuations 
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FIG. 1. Ring with an in-line dot subject to a flux <I> and capacitively coupled to an 
external impedance Z ex f 



are known to reduce the amplitude of the persistent current [7-9] . On the 
other hand, the coherence properties of the ground state have thus far not 
been discussed. It is interesting to ask to what extent zero-point fluctua- 
tions of an electrical environment can influence the ground state of a ring, 
and how they affect the persistent current. This question is also interesting 
in view of the recent discussion on whether or not weak localization prop- 
erties are influenced by zero-point fluctuations [10-14]. This discussion 
focuses on the conductance, a transport coefficient, and not directly on the 
ground state of the system. In this paper, we extend and clarify the results 
presented in [1] . A more intuitive discussion of the results of [1] based on 
a Langcvin approach is given in [15]. Here we provide a detailed deriva- 
tion based on Bethc ansatz results and perturbation theory of the results 
obtained in [1]. In addition, we investigate the capacitance coefficients of 
the ring and the charge-charge correlation spectrum. To be specific, we 
discuss the system shown in Fig. 1, namely a normal ring divided into two 
regions by tunnel barriers and coupled capacitively to an external electrical 
circuit. Each of the two regions can be described by a single potential, and 
the system considered permits thus a simple discussion of its polarizability. 

We turn our attention to the specific case where the external impedance 
Z ext (uj) shows resistive behavior at low frequencies. We consider the model 
discussed in [1], namely an impedance modeled by a transmission line. We 
investigate the persistent current, the charge on the dot and the flux in- 
duced capacitance which is the response of the charge to a variation in the 
magnetic flux through the ring. We identify two different classes of ground 
states as a function of the coupling between the external circuit and the 
ring, namely a "Kondo" type ground state at weak coupling and a symme- 
try broken ground state at strong coupling. In order to understand better 
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the different regimes, we discuss time dependent quantities, in particular 
the power spectrum of the charge fluctuations on the dot. 

The case of a resistive external circuit is of particular importance since 
the ring-dot structure coupled to a resistance exhibits dephasing at fi- 
nite temperatures [15]. In view of the recent controversy concerning the 
question of whether there is dephasing at zero temperature [10, 11, 13] or 
not [12,14], it is interesting to ask what we can learn about the ground state 
of our simple system if it is coupled to a resistive external circuit. The case 
where the dot and the arm of the ring are in resonance has already been 
treated in [1]. Below, we fill in the details left out in the discussion pre- 
sented there and extend it to other thermodynamic quantities as well as to 
the noise spectrum of the charge on the dot. The discussion of the latter is 
of a particular interest since it exhibits decay of the correlations with time 
in the ground state. Here we emphasize an approach which treats the entire 
system in a Hamiltonian fashion. A more intuitive description of this sys- 
tem based on a Langevin equation approach has been given elsewhere [15]. 



2. RESISTIVE EXTERNAL CIRCUIT 

We model the finite impedance Z ext in a Hamiltonian fashion, with the 
help of a transmission line with capacitance Ct and impedance L, see Fig. 2 
and Appendix A. The ohmic resistance generated by the transmission line 
is R = Z ext (uJ = 0) = (L/Ct) 1 ^ 2 . We denote the charge operators on 
the capacitors between the inductances by Q n and introduce the conjugate 
phase operators <j) n fulfilling the commutation relations 

[4>m,Qn] = ie5 mn . (1) 

It is via these charges and phases that the external circuit couples to the 
ring. The charge on the dot is Qd- We define the parallel internal capac- 
itance Ci = Cl + Cr, the external series capacitance C^ 1 = C-f 1 + C^ 1 
and the total capacitances C = d + C e and Cq 1 = C^ 1 + C^ 1 . Then 
the Hamiltonian including all electrical interactions in the circuit reads, 
see Appendix A, 

fy Q 2 d QdQo Qo , fr ( \ 

where Hho is the Hamiltonian of a bath of harmonic oscillators describing 
the transmission line, 
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FIG. 2. Transmission line with impedance Z ex t(t*j) = itoL/2 + \J L/Ct — ui 2 L 2 /4. 



The spectrum of Hho is discussed in some detail in Appendix A. Here, we 
merely state that it is dense. Similar models involving a two-level system 
coupled to a bath of harmonic oscillators with a discrete set of fundamental 
frequencies are investigated in the context of stimulated photon emission. 
A known example is the Jaynes-Cummings model [16]. These models ex- 
hibit recurrence effects such as wave function revival [17, 18]. Because our 
model has a dense spectrum of fundamental frequencies of the harmonic 
oscillators, there are no recurrence phenomena [19]. 

2.1. The total Hamiltonian 

At low enough temperatures, we can describe the ring-dot system in 
terms of just two charge states, namely the state with N electrons and the 
state with iV+1 electrons in the dot. This leads to a two level system [20,21] 
which we discuss briefly in Appendix B, and which is described by the 
Hamiltonian SH^S^ 1 . The total Hamiltonian reads 

ee SHfS- 1 +H C = y *, ^-a x + ^-a z Q + H ho . (4) 

Let us briefly discuss the different terms of Eq. (4) . The Pauli matrix a z is 
related to the charge on the dot by = (e/2)<r z + e(N — N+ + 1/2), where 
eN + is an effective background charge. The detuning he is the energy 
difference that an electron must overcome if it wants to tunnel from the 
arm to the dot. It contains the one-particle energies £d(jv+i) an d e a M ; and 
a term stemming from the Coulomb interactions, Eq. (2). It reads 

= e<j(jv+i) -e aM + ^ (n - N + + ^j . (5) 

The tunneling across the barriers between the dot and the arm is mediated 
by hAo<T x /2, where the tunneling amplitude hA depends on the flux 
$ and is given in Eq. (B.6). The term (e/2Ci)a z Qo is responsible for 
the capacitive coupling between the ring and the external circuit, whose 
Hamiltonian Hho is given in Eq. (3). 

The Hamiltonian H e ^ given in Eq. (4) is the Caldeira-Leggett model 
[19,22]. To make the connection to this model more explicit, we introduce 
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the cut-off frequency uj c = (LC t ) 1//2 of the transmission line and the 
parameter a describing the strength of the coupling to the bath 

(6) 



Rk \ Ci 

where Rk = h/e 2 denotes the quantum of resistance, and R — (L/Ct) 1 ^ 2 
is the zero- frequency impedance of the external circuit. The low energy 
behavior of a system described by the Hamiltonian of Eq. (4) is determined 
by the parameters a, oj c , Ao and e. The limit of a ring which is not coupled 
to an environment is recovered by setting a — 0. This limit corresponds 
to vanishing external capacitances, C\ = C2 = 0, implying that Co = 
0. Before discussing in more detail the different regimes associated with 
different values of those parameters, let us rewrite the coupling between 
the two level system and the bath in a form that is more suitable for our 
purposes. 

We introduce the unitary transformation 

U = cxp(- l -^4> a z ). (7) 



2Q 

With this transformation, we define an equivalent Hamiltonian 

H = UH eff U-\ (8) 

A simple calculation yields H — H TL s+Hho+Hi, where Hho has already 
been defined in Eq. (3), and 

he 

HtLS = -7T<Jz, (9) 



Hj = - 



2 

hA 



' Xi " ' + a ~ cxp \ i lJ~^ 



(10) 



The coupling between the ring and the external circuit is thus mediated by 
the operator exp[— i(Co/Ci)<fio] which increases the charge on the capacitor 
C\ (see Fig. 1) by e Co/Cj. We emphasize that the charge on the capacitor 
C'i is in general not shifted by an integer multiple of the electron charge 
when an electron tunnels between the dot and the arm of the ring. For a 
more thorough discussion, see Appendix C. We investigate the correlator 

P{t) = (ew(i^k{t)}ew(-i^fo{0))) , (H) 



in particular its long time behavior. Here, 4>o{t) = e lHt / h (f>ae lHt / h is the 
time dependent phase operator, and (. . .)o denotes the expectation value 
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taken with respect to the ground state for A = 0, see also Appendix D. It 
is shown in there that the correlator P(t) behaves asymptotically like 

P(T)~(l + »^|) 2tt , (i->oo). (12) 

Eq. (12) is of central importance. It allows to show that at a = 1/2 the 
problem can be mapped onto the exactly solvable Toulouse limit, Sec. 4.2. 
Moreover, it allows for a discussion of the thermodynamic properties of 
the problem using perturbation theory, see Sec. 3.2. Furthermore, the 
correlator P(t) is closely related to the function P(E) encountered in the 
discussion of tunneling through a tunnel barrier shunted by an impedance, 
see [23] and Appendix C. The discussion of P(E) is generally referred to 
as "P(S)-theory" . 

3. GROUND STATE PROPERTIES 

In this section, we discuss ground state properties of the ring-dot struc- 
ture coupled to a resistive external circuit. We are mainly interested in 
the persistent current. However, we also discuss two quantities related to 
the polarizability of the ring-dot structure, namely the charge on the dot, 
Qd and the flux induced capacitance C<[> [20, 21] which is the response of 
the charge on the dot to a variation in the magnetic flux much in the 
same way as the electrochemical capacitance (see, for instance [24]) is the 
response of the charge to a variation in the applied bias. 

3.1. Equivalent models 

The anisotropic Kondo model, see e.g. [25], the resonant level model [26] 
and the inverse square one-dimensional Ising model [27] are equivalent to 
the Caldeira-Leggett model in their low energy properties. For vanishing 
magnetic field, corresponding to our model being at resonance e = 0, those 
models are characterized by two parameters, which may be identified as a 
and Aq/lo c . The half plane spanned by Ao /u) c > and a is divided up into 
three regions by the separatrix equation [1 — a\ = Aq/uj c , where the param- 
eters flow to different fixed points under the action of the renormalization 
group [28], see also Fig. 3. The regime of strong tunneling Aq/uj c > [1 — a\ 
is beyond the scope of this work, and we shall not consider it any further. 
Of the remaining two, the region a > 1, Ao /lu c < a — 1 corresponds to the 
ferromagnetic Kondo model where Aq/lo c renormalizes to zero, whereas the 
region a < 1, Ao /u) c < 1 — a corresponds to the anti- ferromagnetic Kondo 
model that has been solved by Bethe ansatz [29,30]. For large detuning e, 
both cases can be treated by a perturbative expansion in Aq. The details 
of this expansion are treated in Appendix D. 



3.2. Perturbative results 
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a 

FIG. 3. Scaling trajectories in the half plane spanned by a and the normalized 
tunneling amplitude Ao /ui c - If the cut-off is lowered from uj c to lu c — du> c , scaling takes 
place along the lines in the direction of the arrows. The bold lines correspond to the 
scparatrix equation |1 — a\ = Ao/w c , dividing the plane up into the strong tunneling 
(S), the ferromagnetic (F) and the anti-ferromagnetic (AF) Kondo region. This diagram 
has been investigated by Anderson, Yuval, and Hamann [28]. 



At zero temperature and to second order in A , the free energy is F + 5F 
with (see Appendix D) F = —h\e\/2, and we find 



SF 



hu> c 



Ao 

L0 r 



2a-l 



r i 



(13) 



where r(£, x) is the incomplete gamma function. If the detuning e is large 
in magnitude, that is, far away from resonance, Eq. (13) is valid at any 
a up to second order in Ao /w c . In the regime a-1 > Ao/w c (the fer- 
romagnetic regime of the Kondo model, see Fig 3), Eq. (13) is valid even 
for arbitrary detuning e. From Eq. (13), we easily obtain the persistent 
current in perturbation theory 



/ = 



88F He A <9A 



\e\/u 



2 lo c d$ 



r 1 - 2a, 



(14) 



Note that therefore for a > 1, where Eq. (14) is valid for arbitrary detuning, 
the persistent current has a cusp at resonance, see Fig. 4. It is well known 
[31, 32] that the Caldeira-Leggett model exhibits symmetry breaking for 
a > 1, namely (a z ) = dF/d(he) exhibits a finite jump as e crosses zero, see 
Fig. 5. By virtue of Eq. (B.7), (<j z ) is related to the average charge (Qd) on 
the dot. Let us briefly illustrate this relation. For e — > — oo, (<r z ) goes to 
+1, and the charge on the dot approaches therefore an integer multiple of 
the electron charge e{N + 1 — N + ). As e increases and goes towards zero, 
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FIG. 4. The persistent current as a function of the normalized detuning e/ui c for 
a = (solid line), a = 1/2 (dotted line) and a = 1.2 (dashed line). Units are chosen 
such that the persistent current for a = is unity at resonance e = 0. Note the different 
scale for the persistent current at a = 1.2. The parameters are w c = 10 Aq. 



a z approaches as well (unless there is symmetry breaking, in which case 
it approaches a finite positive value). At (cr z ) =0, the charge on the dot is 
a half-integer multiple of the electron charge, (Qd} = e(N +1/2 — N + ), and 
for e — > +oo, the charge on the dot approaches again an integer multiple 
of e. In the following, we discuss for convenience the average (cr z ), bearing 
in mind that the physically more meaningful quantity is the charge on the 
dot (Q d ). 

The average (a z ) is an odd function of e, and for e > it reads in 
perturbation theory 



1 1 / A 



4 V U!, 



101 



E 

n=0 



(-1)' 



(e/w c ) n 



n! (l-2a + n)(2-2a + n) 



•2a- 1 



(15) 



where we have used the series expansion for the incomplete gamma func- 
tion, see Eq. (D.34). In the case of symmetry breaking, a > 1, the width 
of the jump at e — is thus 



£=0- 



e=0+ 



w c y 2(2a- l)(2a-2)' 



(16) 



see also Fig. 5. The jump in (a z ) and the cusp in the persistent current 
have the same origin, thus the cusp is an indicator of symmetry breaking. 
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FIG. 5. The position (a z ) of the electron for a = 1.2 (solid line) showing symmetry 
breaking, compared to the typical behavior of (a z ) at a = (dashed line). The tunneling 
amplitude Aq/u) c = 0.5 is chosen very large in order to make the bending near e = 
visible in the case of symmetry breaking. 



The sensitivity of the charge (Qd) on the dot, see also Eq. (B.7), to 
changes in the magnetic flux $ can be characterized in terms of the flux 
induced capacitance [20,21] 



r = d(Q d ) d(a z ) 
It is an odd function of e as well, and for e > 0, it reads 



(17) 



e A d{A /uj c ) 
2cj c d$ 

- r(l-2a) 



E 

n=0 



(-1)" 



n\ (1 - 2a + n)(2 - 2a + n) 



+ 2a- 1 



(18) 



Like the average electron number, Eq. (16), it exhibits a finite jump at 
resonance 



C*(e = 0+) - Gj,(e = 0-) = e 



A a(A /w c 



1 



u c 9$ (2a-l)(2a-2) 



, (19) 



for a — 1 ^> Ao/uJc- Thus the symmetry breaking can be seen as well in 
the flux induced capacitance, see Fig. 6. Flux sensitive screening in rings 
has recently been observed by Deblock et al. [6] . 
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3.3. Bethe ansatz results 

For < a < 1, we use the Bethe ansatz solution of the resonant level 
model [33]. The low energy properties of the problem depend on three 
energy scales, namely the detuning e, a cut-off D and the "Kondo temper- 
ature" 

T K = A " Q , (20) 

which is the only scale that depends on the magnetic flux <f>. The persistent 
current is calculated from the known expression [33] for (a z ) = dF/d(he) 

, , j r dp M Z+b ^^+^ni/2- l {i- a ) P ) 

{az) ~ 4^/2 p - *0 T(l + iap) ' ( ' 

where z = (e /T^) 2 *- 1- "-' and b = a In a+(l — a) ln(l— a). This Bethe ansatz 
solution is valid for 1 — a > A n /uj c . Now, we make use of the fact that in 
terms of the above-mentioned energy scales, and at zero temperature, the 
free energy reads F = hTxf{s/Tx, Tk/D), where / is a universal function. 
The persistent current I = —c(dF/d&) may thus be expressed in terms of 
(a z ) = dF/d(he) 



<9<f> 



/ dx(a z )(x) -y(a z )(y) 
Jo 



+ I . (22) 



where y = e/T^. The integration constant I is just the persistent current 
at resonance e = and is determined below. The integral in Eq. (22) can 
be performed giving 

hc_ /A \^ dA 

j ~ (_!)« r (\ + n(l - a)) e- nb ^"d-") 



e 2(l-a) 



2(1-^(1-5^)1 
for lnz > —b. Near resonance, more precisely for lnz < —6, it reads 
he 1 / A \ ^ 8A Q 



(23) 



V^2(l-a) \ oj c J 9$ 
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V f + (24) 

t[ n - r (1 + ^-1/2)^) \TkJ 



For a = 1/2, the series in Eq. (23) can be explicitly re-summed 

+ /o- (25) 



1\ He A <9A . 
I [a = - = — In 

2 16 w c 9$ 



I / 16 w c : 



This reflects the fact that the case a = 1/2 corresponds to an exactly 
solvable limit of the Kondo model, the Toulouse limit [34], see also Sec. 4.2. 
It remains to calculate the persistent current at resonance, To- We do so in 
the following paragraph. 

We point out that Tk sets the scale for the transition from resonant to 
perturbative behavior, i.e. for \e\ 3> Tk, the expressions for the persistent 
current from Bethe ansatz and from the perturbation theory must coincide 
up to terms of order A . This observation allows us to determine the cut-off 
D in terms of w c , 

D\ 2a 2T(3/2-a)e- b , , 

- (26) 



y u c J ' ^(1 - 2a)T{\ - 2a)r(l - a) ' 
as well as the integration constant I 

Tic r(l - 2(Tz^y)e" 3 ^ /A \^aA , Tic 1 A 0A 

1 = 77— ^T, 7, x \ -FT I ~~cC7 r 



2 V^(l - a)r(l - 2(T^y) V D J 5$ 2 l-2aw c S$ ' 

(27) 

For a < 1/2, the first term in Eq. (27), behaving like (A /D) a/(1 " a) is 
dominating. For a > 1/2 the second one dominates and behaves as A /D. 
The power law for a > 1/2 is thus the same as one obtains from perturba- 
tion theory at a > 1, using / = —cdF/d$> with F given by Eq. (13). The 
transition from the anti-ferromagnetic to the ferromagnetic Kondo model 
as a crosses 1 is therefore not seen in the persistent current at resonance, 
although it is visible, as stated above, as a cusp at resonance when the 
persistent current is traced as function of the detuning e. 

Leggett et al. [19] investigate coherence properties by discussing the de- 
cay of a state which is prepared such that the electron is e.g. in the dot for 
negative times t < (by applying a large negative detuning e, for example) 
and released at t = (by setting e = for t > 0). They find a decay with 
superimposed oscillations for a < 1/2 and a decay without oscillations for 
a > 1/2. This behavior is interpreted as decohcrcncc setting in at a — 1/2. 
See also our discussion of the charge correlator in Sec. 4. 
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The persistent current, however, remains differcntiablc in e as a passes 
through 1/2 and the ground state retains some coherence even for large a, 
see also Sec. 4. As a function of detuning the cusp at resonance shows up 
only at a > 1. The poles at a = 1/2 in the two terms in Eq. (27) cancel 
each other. Together they give rise to a logarithmic persistent current 

a= ±)*> Hc^ (ln*» - 0.217... ). (28) 

The logarithmic term ln(Ao/w c ) characterizes the transition from power 
law to linear behavior. The poles in the first term in Eq. (27) appearing 
at a = (2n + l)/(2n + 2) for n > are canceled by terms of higher order 
in A . 

The flux induced capacitance C$, see Sec. 3.2, can be obtained directly 
from the Bethe ansatz expression Eq. (21), giving for e > 

e 9A /a$^ (-l)«-i r(±+n(l-a))e-"<> / e y 2 ^ 1 -^ 



e dA a /d® ^ (-l)"- 1 r(i+n(l-a))e-" ft / e \ 
A t[(n-iy. T(l-na) \T K J 



(29) 



for lnz > —b, where z and b have been defined below Eq. (21). For lnz < 
—b, it reads 

e gAp/gg ~ (-1)" T(l + ^)e^ b / £ f "+ 1 
* V5F(l-a) Ao ^ n\ r(l + (n+l/2) T ^)\T K/ 



(30) 



There is no symmetry breaking for a < 1, that is C$ = at resonance, see 
Fig. 6. 



4. CHARGE CORRELATIONS 

In this section, we consider charge-charge correlations on the dot at un- 
equal times. These correlations may serve as an alternative measure of the 
coherence properties of the ground state. We define operator AQ d (t) of 
the charge fluctuations on the dot by 

AQ d (t) = Q d (t)-(Q d ), (31) 

where the time evolution of the operator Q d is given by 



Q d (t)=e l ^ h Q d e- lAt / n . 



(32) 
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FIG. 6. The flux induced capacitance as a function of the detuning e/lo c and for a 
tunneling amplitude Aq/io c = 0.2. The flux induced capacitance is divided by dAo/d& 
such that the graph is independent of the flux <3>, and given in arbitrary units. The solid 
line stands for Cj> at a = 0.3, the long dashed line is Cj> for a = 0.6. The short dashed 
line represents C$ at a = 1.2 blown up by a factor 5 and exhibits symmetry breaking 
in e = 0. 

With Eq. (B.7) the symmetrized charge correlations may be written in 
terms of the Pauli matrix a z 

S QQ {t) ee i({AQ d (t),AQ d (0)}) = j (({* z (t),* z (0)}) -2<^) 2 ) , 

(33) 

where the curly brackets {■, ■} denote the anti-commutator {A, B} = AB + 
BA. In the following we will concentrate on two simple cases, namely a = 
and a = 1/2. 

4.1. Correlations at a — 

The case a = is very easy to treat. In this case the ring is not coupled 
to the external circuit, and it suffices to calculate the charge correlator, 
Eq. (33) for the ground state of the 2 x 2-matrix (he/2)a z — (hA /2)a x , 
see Eq. (4), and we find 

SQQ(t) = T^rk cos (V e2+A o) • ( 34 ) 

The correlations are periodically oscillating with time, and the envelope of 
the oscillations is a constant. The spectral density Sqq(u), defined by 




S qq (cj)= / dte^S QQ (t), 



(35) 
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shows 5-peaks at w = ±\A 2 + A 2 , 



p 2 A 2 
S QQ (u)=Tr e -- A(l 



A 2 



5{u> ^e 2 + A 2 ) + 5(u + y / £ 2 + A§) 



.(36) 



This also reflects the fact that there is no decay of the correlator with time. 

4.2. Correlations at a — 1/2 

In this section, we set % = 1. In the case a = 1/2 the Caldeira-Leggett 
model can be mapped on an exactly solvable limit of the Kondo model, 
the Toulouse limit [34]. Its Hamiltonian describes the hybridization of a 
continuum of spinless electrons with a localized electronic level but without 
electron-electron interactions. This can also be seen from the correlator 
P(t), see Eq. (12), which for a = 1/2 takes the form 

P a =i/2(t) ~ 1 ^ 1 ,,, , (t -> oo). (37) 

l+lOJ c \t\ 

This expression corresponds to the sum over all energies of the Green's 
functions of free electrons in a symmetric wide band with linear dispersion 
and a density of states dn/duj = cj" 1 . We denote the spectrum of the elec- 
trons in the continuum by ui n and their creation and annihilation operators 
by foT and b n , resp. The energy of the localized level is e and its creation 
and annihilation operators are and c. The strength of the hybridization 
is J±. The Toulouse Hamiltonian then reads 



(38) 



Then the Caldeira-Leggett model at a = 1/2 is mapped onto the Toulouse 
limit by identifying the parameters, see [19] 

s <-» e (39) 
A <-» 2~ 1/2 J±. (40) 

Moreover, we identify the Pauli matrix a z in the Caldeira-Leggett model 
with the operator for the occupation of the localized level 



-a x ++$£-\. (41) 



The Green's functions for this Hamiltonian can be calculated in closed 
form. In order to calculate the charge correlator, Eq. (33), we need the 
time ordered Green's function of the localized electron level 

G d (t) = -i (Tc(i)ct(O)) = -i [9{t) (c(i)ct(O)) - 6{-t) (ct(O)c(t))] , (42) 
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where 9{t) is the Heaviside step function. In Fourier space, it takes the 
form 

/OO J 
^e-^G d (t) = [uo-e + iO- E^u;)]" 1 , (43) 
-oo ^ 

where the self-energy £d(w) reads 



S d (w) = -i 



. 7T / J± 



2uj c 



sign a; = — — signw, 



(44) 



where sign a; denotes the sign of u>. We shall see below that V is in fact the 
inverse decay time of the charge correlations. In terms of the parameters 
of the Caldeira-Leggett model, it reads 



2 U>r. ' 



(45) 



Using Eqs. (33) and (41), and with the help of Wick's theorem, wc hnd the 
charge correlations for a = 1/2 to be 



S QQ (t) = e 2 Rc [G d (t)G d (-t)}. 
A rather long but simple calculation yields the spectral density 



with 



l 

47T 

i r 

8ir u 



arctan 



In- 



T/2 



2 2r 

3 ^Tp /(w) ' 



+ arctan 



(46) 



(47) 



e 2 + (r/2) 2 



+ ln- 



r/2 

e 2 + (r/2) 2 



(w - ef + (r/2) 2 (^ + e) 2 + (r/2) 2 
hef us discuss SQQ(t). We have, see also Fig. 7, 



f(u) — > -, for a; 3> max{e, r}, 



(48) 



(49) 



and the Fourier transform of 2r/(e 2 + T 2 ) is simply e~ r '*'. Moreover, 
for t — > 0, the correlation function <Sqq(£) approaches the integral of the 
spectral density 



/•°° 

Sqq(0)=J -^Sqq(uj), 



(50) 
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which, for symmetry reasons, must equal 1/4 for e = 0. A sketch of Sqq(0) 
is given in Fig. 8. The correlation function thus reads for short times 

S QQ (t)~S QQ (0)e- r ^, (|t|«min{ £ - 1 ,r- 1 }). (51) 

For long times, we use the stationary phase expansion and get 

S QQ (t)~~ — (r/2)2 A ' 2 , (i-Kx>). (52) 

QQK ' 7T [ £ 2 + (r/2) 2 ] 2 v ' y ' 

Thus the correlations decay exponentially with rate T at small times and 
algebraically for large times. Moreover, since Sqq(uj = 0) = 0, the integral 
of Sqq (t) over all times must vanish, therefore Sqq (t) takes negative values 
at intermediate times, see Fig. 9. For e ^ the correlations are still 
oscillatory in time. 

We have discussed here only the cases a = and a = 1/2 which corre- 
spond to no coupling at all between the ring and the electric circuit (a = 0), 
and a relatively strong coupling (a = 1/2). At a = 0, the charge correla- 
tion spectrum reveals simply the coherent charge oscillations between the 
dot and the arm with a sharp frequency. At a = 1/2, charge transfer still 
occurs, but it is now entirely characterized by the (flux-sensitive) relaxation 
rate F and the detuning. Instead of the periodic charge transfer which ex- 
ists for a = 0, the charge transfer has a much more stochastic character. 
Together with the results for the persistent current and the polarizability of 
the ring obtained in the previous section, the charge correlation spectrum, 
discussed here, provides an additional measure of the coherence properties 
of the ground state. 



5. CONCLUSIONS 

We have investigated a simple model of a mesoscopic normal ring cou- 
pled to an external electrical circuit. We have shown that the zero-point 
fluctuations of the external circuit influence measurable quantities like the 
persistent current /, the flux induced capacitance C$ , and the polarization 
of the ring. The reduction of the persistent current as well as the onset of 
symmetry breaking, observed in the polarization and in the flux induced 
capacitance indicate that the coherence of the wave function in the ring 
is reduced by fluctuations in the external circuit. While such a suppres- 
sion of quantum coherence is not too surprising at finite temperatures, we 
have demonstrated in the present paper that it persists even in the extreme 
quantum limit of zero temperature. 



APPENDIX A 
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CO 

FIG. 7. The function f(ui) describing the deviation of the spectral density Sqq(lu) 
at a = 1/2 from the spectral density of an exponentially decaying correlator. The 
frequency ui is measured in units of T. 



0.3 




1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

1 2 3 4 5 

e 

FIG. 8. Instantaneous charge correlations as a function of e, which is measured 
in units of T at a = 1/2. The function is symmetric under e — ► — e. The correlations 
SQQ(t = 0) (denoted by 5(0) in the figure) are measured in units of the squared electron 
charge e 2 . 
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0.25 



e = 

e = r 

e = 5r 



exp[-n]/4 







1 



2 



3 



FIG. 9. Time dependence of the charge correlations Sqq(*) (denoted by S(t) in the 
figure) for different values of e in units of the squared electron charge e 2 at a = 1/2. 
The line showing the exponential decay 1/4 cxp(— Tt) is given for comparison. For e ^ 
the correlator exhibits damped oscillations. 



In order to understand the dynamics of the bath of harmonic oscilla- 
tors describing the transmission line shown in Fig. 2, we have to discuss 
its properties in some detail. Its impedance is of the form Z ext {uj) = 
ituL/2 + (L/C t - uj 2 L 2 /4) 1/2 , see e.g. [35]. This appendix is devoted to 
the discussion of the Hamiltonian of the transmission line and its eigen- 
states. We denote the charges on the capacitors between the inductances 
by Q n and the potentials by V n , n = 0, 1, 2, . . .. The charge Qo plays a 
special role, see Fig. 1, because it is the charge sitting on the capacitor Ci; 
Vo is the corresponding potential. Similarly, Qoc and Vac are the charge 
and the potential on the capacitor C 2 • It is via these charges and potentials 
that the external circuit couples to the ring. Furthermore, we introduce 
the phases which satisfy the equations 



The canonically conjugate variables to the charges Q n are not the phases 
4> n but rather the generalized fluxes (fi/e)<f> n . We choose the dimension- 
less quantities 4>„ for later convenience. In terms of these variables the 
Hamiltonian of the transmission line reads 



The Hamiltonian of the transmission line 



= £ (v - v ) 



(A.l) 




h\ (4> n +i - 

e 2L 




(A.2) 
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which is canonically quantized by replacing the charges and phases by op- 
erators Q n and <j) n and by imposing the commutation relations \<f> m , Q n ] = 
ie 5 mn , giving Eq. (3). We note that Q n plays the role of a momentum, and 
4> n plays the role of a position. In particular, the quadratic form involv- 
ing the charge Q n is positive definite. Due to a theorem of linear algebra, 
any pair of quadratic forms, one of which is positive definite, can be si- 
multaneously diagonalized. Therefore, the Hamiltonian Hho of Eq. (3) 
can be diagonalized by a unitary transformation mapping Q n onto Q(x), 
and <f) n onto 4>{x). The transformed operators Q{x) and 4>{x) satisfy the 
commutation relations [<j>{x),Q{y)\ = ie5(x — y). In this new basis Hho 
reads 




Hao = h dx \W + \-e) r in ItJ • ^ 

In terms of the operators Q(x), the charge on the dot Q is given by 

Qo = V2j dx cos (^) Q{x), (A.4) 

and analogously, we have 

$o = V2 J dx cos {^^4>{x). (A.5) 

Eq. (A.5) will be useful in Appendix D in order to calculate the correlator 
-P(t). The diagonalized Hamiltonian Eq. (A. 3), being a sum over harmonic 
oscillators, can be written in terms of the bosonic creation and annihilation 
operators aj. and d x which are related to the charges Q(x) and the phases 
4>(x) by 

fit) , (A.6) 

2 (ox + 4). ( A -7) 

with \ x = h(CtfL) 1 ^ 2 sm(Trx/2), and which obey the commutation rela- 
tions [a x , a) y ] = 6(x — y). In terms of the operators a x and a x , the Hamil- 
tonian of the transmission line, Eq. (A. 3) reads 




Hho = J dxhu x [a) x d x + , (Ai 
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the cigenfrequencies u> x being 




(A.9) 



We note that the spectrum Eq. (A.9) is linear at low frequencies with a 
density of states dx/dcu = (2/tt)u>~ 1 , where lo c = (LCt)^ 1 ^ 2 , and that it is 
cut off at a frequency ui c . We shall need Eqs. (A. 8) and (A.9) in Appendix D 
where we write the partition function of the Caldeira-Leggett model as a 
path integral. 



APPENDIX B 

Degrees of freedom in the ring 

The Hamiltonian He of the Coulomb interactions, Eq. (2), has to be 
augmented by a Hamiltonian describing the electronic degrees of freedom 
inside the ring. We are interested in the effect of the long range nature 
of the Coulomb interactions, therefore we disregard charge redistribution 
and correlation effects and assume the potentials in the dot and the arm of 
the ring to be constant over the extent of the respective subsystem. As a 
consequence, the electronic wave functions arc independent of the number 
of electrons in the dot and the arm. Moreover, we consider spinless electrons 
here. It has been demonstrated in [20,21] that the inclusion of spin docs 
not change the calculations a great deal although it somewhat changes the 
physics of the problem. In particular, the Kondo effect for an isolated ring 
sets in only for a tunneling amplitude hA that is comparable to or larger 
than the mean level spacing in the dot or the arm, whereas in this paper 
we assume that UA is much smaller than the mean level spacings. For 
a discussion of the Kondo effect in normal metal rings, see the articles of 
Eckle et al. [36] and of Kang and Shin [37] . We denote the creation operator 
for an electron of energy edj in the dot by dj and the creation operator for 

an electron of energy e a k in the arm by a\. The tunneling through the 
left barrier is described by a hopping amplitude tj k , the tunneling through 
the right barrier by t^ k . The total tunneling amplitude includes the phase 
picked up by an electron whose wave function encloses the magnetic flux 
$ and reads 

t jfc ($) = tf k + tf k exp (2™J^ . (B.l) 

where $ = hc/e is the elementary flux quantum. As the dependence of the 
Hamiltonian on the magnetic field is shifted into the boundary conditions 
by Eq. (B.l), the tunneling amplitudes tj k and t^ k can be taken to be real. 
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Their relative sign depends on the number of nodes of the wavcfunctions 
they couple. It is positive if the sum of the number of the nodes of the 
wavefunction in the dot with energy and of the number of nodes of the 
wavcfunction in the arm with energy e a k is even, and negative otherwise. 
This makes up for the parity effect in the persistent current. Now, we can 
write a noninteracting Hamiltonian for the electronic degrees of freedom 



E 



jk 



[t ifc (*)dt 



a k + h.c. 



(B.2) 



All the many-body interactions are taken into account by He of Eq. (2). 
The operator for the charge on the dot becomes 



Qd = e^d]d - eN + 



3 



(B.3) 



with eN + being an effective background charge on the dot. Next, we make 
a crucial approximation, which allows us to write H e and Qd as 2 x 2- 
matrices. 

We shall assume in the following that the tunneling amplitudes through 
the left and the right barriers are much smaller in magnitude than the 
level spacing in the dot and the level spacing in the arm. We note that the 
number of charge carriers in the ring is conserved. Following Stafford and 
one of the authors [20,21], we consider hybridization between the topmost 
occupied electron level in the arm and the lowest unoccupied electron level 
in the dot, e a M and e^N+i) only. To simplify the notation, we denote 
the tunneling amplitudes between the levels e a M and e^jv+i) by th for 
tunneling through the left barrier and by t R for tunneling through the 
right one, and we assume tL and tn to be positive. The Hamiltonian 
for the electronic degrees of freedom can be written in terms of the Pauli 
matrices a z and <j± = (<r x ± iv y )/2 



jjeff _ e d(N+l) - CgM 



t L ± t R exp 



c 2 - 



t L ± tji exp ( 2iri 



-2m — 
$ 



En 



(B.4) 



where E = Ylj=i e dj + Y<k=i e "fe + ( e "M + e<i(JV+i))/2. The terms in a± 
describe the tunneling between the states £d(N+i) an d e a M- The relative 
sign of tL and tR is positive if the number N + M of electrons in the ring 
is odd, and negative if it is even, as discussed in the previous paragraph. 
The unitary transformation H^p SH^pS~ x where S — exp(—i<j z 'd/2) 
and tantf = -t R sin(27r$/$ )/(^L + t R cos(27r$/$ )), leaves the term in 
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a z unchanged and transforms the tunneling term as follows 



S { \t L ± t R e 2 ^^" 



?jA 



(B.5) 



with 



= Jt 2 L +t 2 R ±2t L t R cos [ 2tt— 



(B.6) 



the positive sign applying again for an odd number of electrons, the negative 
one for an even number of electrons in the structure. The charge on the 
dot, Qd, reads in the two level approximation 



Q d =-<J z +e[N - N+ + \) . 



(B.7) 



APPENDIX C 

Relation to P(i?)-theory 

In this section we discuss the relation of the results obtained above to 
the "P(i?)-theory" [23]. Usually, the P(_E)-theory is employed in order 
to calculate the tunneling rate through a barrier which is shunted by an 
impedance Z s h U nt{w) in the Golden Rule approximation. The system dis- 
cussed here is somewhat special in the sense that the impedance Z ex t(u>) 
does not act as a shunt, but is coupled only capacitively to the tunnel bar- 
rier, see Fig. 1. Therefore, the electrons in the ring cannot flow through the 
impedance as is the case for a shunt impedance. We can see the difference 
to ordinary P(_E)-theory more clearly by considering the Hamiltonian Hj, 
Eq. (10), responsible for the coupling between the ring and the transmis- 
sion line. The operator exp[— i(Co/Ci)<j>o] appearing in Eq. (10) is a charge 
shift operator. It differs from the charge shift operator that is usually en- 
countered in the context of tunneling in the P(i?)-theory in that it does 
not change the charge on the capacitor C\ by an integer multiple of the 
electron charge but by a fraction (Co/Ci)e, 

e-HCo/c^oQ^Co/c^o =q q _ Co e . (CJ) 

c« 

The difference is due to the fact that the standard P(i^)-theory does not 
consider <f> but an operator <f> which is canonically conjugate to the charge 
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Qd on the dot. This is no problem in an open system, where the conjugate 
operator <j> indeed exists. In our closed system with only two charge states, 
however, there is no conjugate operator to Qd, therefore we stick to the 
operator exp[i(Co/Ci)<po\- The coupling to the external circuit modifies 
the tunneling such that the persistent current is reduced, as shown in the 
present paper. The modification is described by the correlator 

K{t) ^ (§) (^(*)^(°)-^(°)) > ( C - 2 ) 

where (. . .)o denotes an expectation value with respect to the ground state 
for A = 0, see also Appendix D. The pre-factor (Co/Cj) 2 is character- 
istic for our problem and is due to the fact that we are considering the 
phase operator (j) of the external circuit. The correlator K{t) is related to 
the function Pit) given in Eq. (11) by Pit) = exp[.ff(f)], and the Fourier 
transform P(E) of P(t) has given the theory its name. For the purpose of 
calculations it is more convenient, however, to stay with the correlator K(t). 
It has been mentioned above, see Eq. (12), that K(t) <~ — 2aln(l + iw c |t|) 
in the long time limit, where w c is the cut-off frequency, and a is identical 
to the coupling parameter given in Eq. (6), which determines the suppres- 
sion of the persistent current, Eq. (27). In P(.E)-theory, a general rule for 
calculating K(t) at finite temperature is given, 

In the zero temperature limit, this formula yields the same long time behav- 
ior as the one obtained above for a special realization of the transmission 
line. 



APPENDIX D 
Thermodynamic properties of the Caldeira-Leggett model 

We want to write the Caldeira-Leggett Hamiltonian H = Htls + Hho + 
Hi in a form that is suitable for the path-integral formalism. For the bath 
of harmonic oscillators, described by Hho, this has been done in Eq. (A. 8). 
The Hamiltonians containing the Pauli matrices a z and a± are translated 
into the operator formalism by introducing the fcrmion operators c' , c, 
creating and annihilating an electron the dot, and and c, creating and 
annihilating an electron in the arm 

H' TLS = a±($c-$$), (D.l) 
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H i = T~ 



c^cexp ( i^-(f>o J + c^cexp [ — i— ■ 



• (D.2) 



For the phase operator (j> , responsible for the coupling between the ring 
and the external circuit, we have 

- e f 1 , cos(7rx/2) , „ + N , . 

ifa = -j= I dx ±=L ' (a x - fit) . (D.3) 



t Jo 

Thus the Hamiltonians H' TLS , H\ and i?^o can ah be expressed in terms 
of the boson operators d x and a x , and in terms of the fermion operators c, 

ct and c, c 1 ^ . 

In order to identify the Hamiltonians H and iJ' = H' TLS + i/^fo + ^7 wc 
have to exclude the case where both the dot and the arm level are occupied 
and the case where both are empty, in other words, we have to exclude the 
subspaces where the operator h = c) c + etc takes the eigenvalues and 
2, resp. However, the Hamiltonian H 1 docs not change the number of 
electrons in the ring, 

[H',h]=0. (DA) 

Therefore, the partition function Z' = Tr exp(— f3H') is of the form Z' = 
Zq + Z\ + Z2, where the subscript denotes the eigenvalue of h, with Z\ 
being equal to the partition function Z = Tr exp(—0H) of the Caldeira- 
Leggett model. Moreover, on the subspaces with n = and n = 2 the 
two level system and the bath of harmonic oscillators do not interact, and 
on the same subspaces H' TLS = 0. Thus, Z n = Z 2 = Zho with Zho = 
Tr ex.p(~/3HHo), whereas Z\ — ZhoZtlsZi- Here, Ztls is the partition 
function of H' TLS restricted to the subspace where n = 1, or cquivalcntly 
the partition function of Htls, which reads Ztls = 2cosh((3he/2). The 
partition function Z' reads thus 

Z' = 2Z HO (l + Z.cosh^ , (D.5) 

and the partition function Z of the Caldeira-Leggett model becomes con- 
sequently 

Z = 2Z ffo Z/cosh tL— . (D.6) 

Eq. (D.6) gives us the relation between the Caldeira-Leggett model and the 
Hamiltonian Eqs. (A. 8) and (D.l, D.2). We shall exploit this relation to 
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evaluate the partition function of the Caldeira-Leggett model in the path 
integral formalism. 

For this purpose, we determine the partition function Z 1 — Tr exp(— /3H') 
in the path integral formalism. The time dependent creation and annihi- 
lation operators of the harmonic oscillators become complex variables 



a x (t) = e iflt ' h a x e 



iHt/h 



e iHt/Htf e -iHt/n 



a x (t), 



(D.7) 
(D.8) 



and the operators for the electronic degrees of freedom become Grassmann 
variables 



c(t)^c(t), t{t)^c(t), 

C + (t) !-» C*(t), j t (t)HC*({). 



(D.9) 
(D.10) 



We write the imaginary time (r = it) action functional S = Stls + Sho + 
Si with 



Stls 
Sho 
Si 



dr 



dr 



c*(t) 



d_ 



c(r) 



dxa*{r) 



Ao f 
2 J 



dr 



Co 



c*(t)c(t) exp ( i— 0o (t) 



.Co 



+ c*(t)c(t)cxp ( -i — o (r) 



(P.ll) 
(D.12) 

(D.13) 



With these definitions, the partition function can be written as a path 
integral [38] 



VcVc* VcVc* VaVa*e- s , 



(D.14) 



where Va and 2?a* stand for integration over all complex variables cl x (t), 
a*(r), and T>c, . . . etc. stand for integration over the Grassmann variables 
c(t), . . . etc. The expectation value of a function A(c, c* , c, c* , {a x }, {a*}) 
for vanishing tunneling amplitude is given by 



VcVc* VcVc* VaVa* 



1 

= Z~o 

x A(c, c*,c,c*,{a x }, {a* x })e-^ s+SHO \ 



(D.15) 
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where Zq denotes the partition function for vanishing tunneling amplitude 



J 



VcVc* VcVc* VaVa* e -( s TLs+s HO ) 



thus the partition function Zj reads 



= (^)o- 



(D.16) 



(D.17) 



This expression can be written as a series expansion in even powers of 
Ao. It has been shown in [15] that the partition function thus obtained 
is formally equivalent to the partition function of the anisotropic Kondo 
model [28] . Below, we calculate the free energy to second order in Ao from 
Eq. (D.17). This perturbative result is valid for large detuning e only, if 
a < 1, and even at resonance s = in the case of symmetry breaking, 
a > 1. 

Before calculating the free energy, we introduce the Green's functions 
needed. The electronic degrees of freedom are described by the standard 
Matsubara Green's functions G(r — r') = (c(t)c* (t')) and G(t — r') = 
(c(t)c*(t'))o, which are anti-periodic with period (3, that is, G(r + /?) = 
-G(t) and G(t + (3) = -G(t). On the domain -(3 < r < (3 they read [39] 



G(t) = e- eT ' 2 



G(r) = -e £T/2 



6(r)- 
9(-r) 



1 



e 0e/2 + 1 
1 



e f3e/2 + l 



(D.18) 

(D.19) 



where 0{t) is the Heaviside step function. The relevant operator <j) of the 
bath, see Eq. (D.2), has its dynamics described by the correlator 



e i(C o /C i )0o(r) e - 4 (Co/C I )0o(O) 



(D.20) 



where K{t) = (G /G 4 ) 2 (M T )M^) - <£o(°))o- This correlator, exp[K(r)] 
is equivalent to the correlator P(t), Eqs. (11, 12), evaluated at imaginary 
times. The correlator K(t) is easily computed for the transmission line, 
Eq. (A. 8). On the domain — /3 < r < (3 the Green's functions of the bosons 
a x read [39] 



(a x (T)a* y (0)) =S(x-y) ( 



6{t) 



CPU* 



1 



(D.21) 



Being of bosonic nature, they are periodic in r with period (3. At zero 
temperature, (3 — > oo, we may drop the second term in square brackets in 
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Eq. (D.21) and obtain the correlator for <f>o(t) 

(Mr)M0)) o =clj o dx huJx ■ (D.22) 

We linearize the spectrum of the harmonic oscillators 

huj x — * —hid c x, (D.23) 

see also the discussion of the spectrum of the harmonic oscillators following 
Eq. (A. 9). The cosine in the integral in Eq. (D.22) is replaced by 1, and 
the integral is cut off at a frequency u c . We find 

<<Mt)<M0)) o ~ — / — e-Me-'M (r ^ oo), (D.24) 

tiK Jo w 

where we have used the quantum of resistance Rk = h/e 2 and the resis- 
tance of the transmission line R = (L/Ct) 1 ^ 2 - This integral is infrared di- 
vergent, a disease remedied by subtracting the correlation function (0§(O)). 
We obtain the correlator K(t) 

K(r)~-2a — 1 - e^ 7 "' e"^ = -2a In (1 + w c |t|) (t-»oo), 
Jo u L J 

(D.25) 

thereby defining the dimensionless parameter a describing the strength of 
the coupling between the two level system and the bath 

R ( Co x 



From K(t), we can calculate thermodynamical properties from the parti- 
tion function. We shall do so in the next paragraph. 

From Eq. (D.6), we can calculate the free energy by the relation — (3F = 
\vlZ. In fact, the partition function of the bath, Zho, only gives rise to 
an additive constant which wc shall neglect below. We therefore write 
F = F + SF with 

F a = -k B T\nZ TLS ^-^- (T^O), (D.27) 

where k B is the Boltzmann constant, and T is the temperature. The cor- 
rection to the free energy due to the presence of the bath is 



SF = -k B T\nZi. 



(D.28) 
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We shall use the cumulant expansion of the expression for Z/, Eq. (D.17), 
to calculate SF to second order in A 



Zi = {e- S ') = cxp [- (Sj) + \ ((S]) Q (Stf) ... 



. (D.29) 



As (Si) a = 0, we only need (Sj) . We find 



(S]) = -A 2 J^ dr cIt'G(t-t')G(t' -T)e K(T '- T \ (D.30) 

Now, we use Eqs. (D.29, D.28), to calculate the correction SF to the free 
energy to second order in Ao- We assume for convenience that e > 0. The 
case e < follows immediately from the symmetry of the free energy 



F(-e)=F(s). 



(D.31) 



For e > 0, we insert the zero temperature limits of the fermion Green's 
functions Eqs. (D.18, D.19) into (Sj) , Eq. (D.30), and obtain 

and take the limit of zero temperature, (3 — ► oo. Taking into account the 
symmetry relation, Eq. (D.31), this gives 
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(D.33) 



where r(£, a;) denotes the incomplete gamma function. Using the series 
expansion for r(£, x) 



T(C,x)=T(C)-J2 



(-l) n a^ +n 
n!(C + n) 



(D.34) 



n=0 

and the shorthand notation x = e/w c this becomes 

2 



5F = _^ £ ^Ao ] eW 



1 2a 



- 1 r(l-2a,H)-^-^ 



(-i)"N r 



n=0 



n!(l - 2a + n) 



This expression converges at resonance, e = 0, for any a > 1/2 



5 F(e = 0) 



4(2a- 1) 



Ao 



(D.35) 
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This result coincides, for a > 1, with the numerical value obtained by 
integrating the renormalization flow equations for the ferromagnetic aniso- 
tropic Kondo model [28]. The expectation value (a z ) = dF/d(Tie), how- 
ever, diverges for a = 1 in perturbation theory. Moreover, in perturbation 
theory, the condition —1/2 < (<j z ) < 1/2 is fulfilled only for a > 1 and 
Ao /u) c -C a — 1, which corresponds to the ferromagnetic regime of the 
anisotropic Kondo model, see Fig. 3. It is in this domain we believe the 
perturbation theory to apply for all e. 
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